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some recently proposed entangling projections. 
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Introduction 

After the pioneering works in 1974 and 1976 by 
Davies ([U, a huge amount of physical information 
about the irreversibility and the evolution of open quan- 
tum mechanical systems has been gained. Alicky [3] 
showed in 1977 that these efforts were deeply con- 
nected to the celebrated "Fermi's Golden Rule", that 
now had become mathematically consistent. The con- 
ceptual importance of these works is clearly not only 
academic, as the need for a better understanding of ir- 
reversible processes has never been more urgent. Today, 
so many nanotechnologies are pushing devices towards 
limits where neither quantum phase coherence, nor dis- 
sipation/dephasing, can be neglected 0, H, [l]. Many at- 
tempts to improve the theory have been made since then 
(see for example 0, 0); but despite the compelling need, 
no substantial, fundamental progress, directly applicable 
to nowadays technologies, has been made so far. 

To be more specific, the problem is to understand the 
dynamics of a subsystem of interest, when the global sys- 
tem is fully coherent. In many cases, the dynamics of 
the global system can be splitted into a part that leaves 
the subsystem invariant, plus an interaction between the 
subsystem and the remaining "uninteresting" degrees of 
freedom. The problem then arises whether or not the 
subsystem can be given a markovian, possibly dissipa- 
tive, dynamics as a consequence of its interaction with 
those degrees of freedom. This is generally impossible 
of course, but it has been shown since the '70s that a 
positive answer can be given when the amount of unin- 
teresting degrees of freedom is huge, and the interaction 
is made small. This last condition is referred to as the 
"weak-coupling limit". 

In [l[ the author was able to give a solid physical 
model of a discrete " atom" (system A) interacting with a 
fermionic particle reservoir at thermal equilibrium (envi- 
ronment B). In that case, the subsystem was made of un- 
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entangled pairs of atomic states (also referred to as "den- 
sity matrices" ) and a fermionic thermal equilibrium state, 
while all the entangled pairs just constituted the remain- 
ing, uninteresting, degrees of freedom. The model was 
of high conceptual importance, as the markovian (and 
dissipative) dynamics for the subsystem was shown to 
guarantee the state positivity at all times. This fact gave 
just enough internal robustness to the model as to be of 
invaluable practical use, for at any time, the subsystem 
evolution could be given a strong physical meaning. But 
the system A had to be finite dimensional, or at least, 
its unperturbed hamiltonian was forced to have discrete 
spectrum. 

Since then, much research has been (and currently is) 
focusing on the so called "projection techniques". The 
name "projection" comes indeed from the fact that, as 
explained brilliantly in a certain projection superop- 
erator Pq is introduced: Pq "is the mathematical expres- 
sion for the idea of the elimination of degrees of freedom 
from the complete description of the states of the total 
system" . The effort is here to understand what could be 
the basic requirement for choosing the subsystem/Po in 
order to guarantee a positive subsystem dynamics. Much 
information on the general properties of some bipartite 
subsystems has thus been gained d, [13, UII , which is rel- 
evant in the field of (possibly entangled) Open Quantum 
Systems, where one studies the dynamics of a "small" 
system A interacting with a "big" environment B (not 
necessarily at thermal equilibrium). Usually, one sup- 
poses to be given a generic " phenomenological" marko- 
vian and dissipative dynamics for the global system, and 
writes down a (markovian dissipative) master equation 
for the subsystem. But, as we shall see in the following, 
the strong interplay between the projection (i.e. subsys- 
tem) and the unperturbed global dynamics, severely re- 
stricts the possibilities for Pq, or on the other side, for the 
global system dynamics. For this reason, although these 
models are much relevant to nowadays technologies, they 
all still suffer from the original restriction, that system 
A should be finite dimensional, or at least it should have 
discrete spectral properties, just like the original Davies' 
model [l|. We have to say here that alternative impor- 
tant methods have been addressed since a long time (see 
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for example [12] )i and still are (see [Tl|), which are based 
on stochastic techniques: these models are of course of 
enormous practical use, but since some randomness has 
to be put by hand, they are somehow of less fundamental 
nature. 

Unfortunately, the possibility to obtain a "completely 
positive markovian subsystem dynamics" — also referred 
to as Quantum Dynamical Semigroups, or QDSs — 
through the weak-coupling limit procedure, becomes no 
more available when the system is infinite dimensional. 
In that case indeed Davies showed that a markovian ap- 
proximation, for the subsystem exact projected dynam- 
ics, could still be achieved [1], but positivity could not 
be shown in general, if not up to finite times: not even 
for the old and "safe" partial tracing over the thermal 
bath's degrees of freedom. Although more general, the 
theory did not share the enormous success of the previous 
one among physicists, precisely because of this serious 
limitation. For example, all of the steady state analysis 
became, physically, completely meaningless. 

In this work we follow the path set by Davies, in sup- 
posing the global system is undergoing a fully coherent 
(hamiltonian) dynamics, and propose a new master equa- 
tion for the subsystem, on the basis of time-symmetry 
considerations. At this point, we will introduce a new ob- 
ject, called "completed collision time", which represents 
a safe overestimation of the subsystem relaxation times, 
as well as a safe underestimation of the global system 
recurrence time. This will be a key physical step, that is 
amply justified from the following well known fact (see 
for example [l^, [ij, among many others) : the " semi- 
group" or markovian approximation is physically mean- 
ingful only when neither we observe the subsystem at 
short times (where probabilities of excited states have 
a parabolic time dependence), nor we observe it long 
enough to see Poincare recurrence cycles. Indeed these 
cycles are present, and explicitly appear in the subsys- 
tem exact dynamics, whenever the global dynamics is 
hamiltonian. This time will scale with the coupling con- 
stant, and will serve us to "dynamically diagonalize" our 
time-symmetrized master equation, thereby eliminating 
in symmetric fashion most of the off-diagonal terms in the 
generator. We will then be able to show that our mas- 
ter equation will correctly describe the exact subsystem 
dynamics in the weak-coupling limit, it will be defined ir- 
respective of the subsystem dimensions or spectral prop- 
erties, and it will guarantee positivity for the dynamics, 
for a plethora of different projections/subsystems. The 
number of possible projections is indeed huge enough 
to incorporate the partial trace over the bath, together 
with some more recently proposed "entangling" projec- 
tions and some others by us, that furnish a new 
Quantum version for the Fermi's Golden Rule (QFGR). 

This deserves a comment: there are many proposed 
quanturn generalization of the well known Fermi's Golden 
Rule 0, [3 1 which are robust and physically and mathe- 
matically meaningful. All these generalizations consider 
a bipartite system (this in fact is the only way to pro- 



pose the standard Fermi's Golden Rule in mathematical 
terms, see [1]): nevertheless, the original idea by Fermi, 
although mathematically ill-posed, did not come from 
any sort of (rather sophisticated) projection on some bi- 
partite system. Stated in modern terms, Fermi's idea was 
rather to take a global system and project on the space 
of density matrices, that are diagonal in the basis of the 
unperturbed hamiltonian. His motivations referred to 
the fact that the system eigenvalues are "robust" against 
dissipation, and thus constitute the relevant degrees of 
freedom. Just along these lines, no environment will be 
present in our model of the QFGR, and the relevant de- 
grees of freedom will be " box-diagonal" density matrices 
for a " box-diagonal" unperturbed hamiltonian. The pos- 
sibility to obtain a semigroup (a dissipative and marko- 
vian dynamics) comes exactly by the fact that the unper- 
turbed hamiltonian has continuous spectrum, thus con- 
ferring our QFGR version an autonomous relevance. 

It is worth noting that the explicit form for our master 
equation in this case will share some very similar and in- 
teresting formal properties with some recently appeared 
master equations comingfrom the theory of Entangled 
Open Quantum Systems [1, [13] ■ These, as ours, are hnear 
homogeneous equations for a " composed" density matrix, 
which turn into a linear homogeneous system of equa- 
tions for each component. Needless to say, these kind of 
equations are currently showing an enormous potential, 
and will be shown here to be taken into account by our 
general theory. But contrary to our QFGR, these models 
are all based on bipartite systems. 

Two final introductory remarks. First, we have to say 
that, contrary to what is sometimes believed, a "diago- 
nal" superoperator does not imply that the dynamics is 
free from quantum effects: first of all, "diagonal" refers 
to some basis, say for example that induced by some un- 
perturbed hamiltonian Hq. Then, talking in Heisenberg 
picture, it is true that a diagonal generator will leave a 
subset of observables — those commuting with Hq — in- 
variant (and this would be a remarkable property!). But 
the subsystem needs not be formed only by such a subset, 
and nontrivial commutation relations could well be ob- 
served within the subsystem evolution, even for "diago- 
nal" generators. With respect to this, our theory is highly 
non-trivial, as it fully keeps track of quantum effects and 
does not, to any extent, constitute a sort of semiclassi- 
cal limit. This is at variance, for example, with many 
widely used Green functions truncation schemes, where 
one takes "n-point" quantum correlations functions into 
account, but only up to some finite n (from this point of 
view, the weak-coupling and the Green function theories 
are applied in completely different contests, a fact that is 
often misunderstood). Rather, the fact that a generator 
is diagonal in the unperturbed basis will be shown to be 
strictly linked with its algebraic properties, like that of 
the positivity of its generated dynamics. 

Finally, we have to say that we have made explicit use 
of Hilbert spaces instead of the (much more natural en- 
vironment of) Banach spaces. This is justified by the 
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fact that the Banach space TiTi) of trace class opera- 
tors on a Hilbert space Ti can be equipped with a scalar 
product, and the resulting completion is a Hilbert space 
(that of of Hilbcrt-Schmidt operators on Ti). Our choice 
only conies from simplicity in the proofs involved, as we 
can make use a huge amount of spectral theory, but of 
course is still a rather serious limitation. Nevertheless 
we feel that our message could be extended to Banach 
spaces, and we are currently obtaining promising results 
in this direction. This indicates that the physical content 
of what we discuss here should survive the mathematical 
formalism, substantially unmodified. We also want to as- 
sure the less mathematically inclined (or more physically 
motivated) reader, that our theorems, although mathe- 
matical statements, will be critically discussed within a 
fully self-contained physical contest, that will offer, so 
to say, the physical arguments that led us to the proofs. 
These will instead be reported in the appendices, and are 
directed to the more mathematically oriented readers. 

The paper is organized as follows: in the next Section 
we begin with a warm up to clarify notations and state 
some physical example to motivate our theory, which we 
then present. In Section [IT] we discuss a new time sym- 
metrized markovian generator and make contact with 
some preliminary results by one of us, reported in p^ : 
m Section [ml we propose our dynamically time averaged 
generator; finally in section HVl we find some fairly gen- 
eral class of projections that make our generator into a 
Lindblad form [13], thus obtaining a generator of a (com- 
pletely positive markovian) Quantum Dynamical Semi- 
group. We will also discuss some (already cited) impor- 
tant applications, such as our Quantum version of the 
Fermi's Golden Rule and some recently appeared entan- 
gling projection [9]. Then we shall conclude, and report 
our proofs in the Appendices. 



I. GENERAL FRAMEWORK 

Let us briefly run through some well known examples, 
that will serve us to better state what we aim at, clarify 
notation and motivate the general theory. 



A. Partial Trace 

This is also referred to as "tracing away the bath de- 
grees of freedom" : let H = Ha (^Hb the tensor prod- 
uct of two Hilbert spaces, representing a coupled system- 
environment pair. Let B be the space of trace-class oper- 
ators on Ti. and Bq the space of trace-class operators on 
Ha- These are Banach spaces, but they can be completed 
with the scalar product {A, B) = Ti{A'^ B) to the Hilbert 
spaces (respectively H.S.{H) and H.S.{Ha)) of Hilbert- 
Schmidt operators on the respective Hilbert spaces: we 
shall always work within these spaces, which we shall de- 
note with the same letters B and Bq. Now define the lin- 
ear map, called the "partial trace" , Pq : B ^ Bq uniquely 



determined by 

Tr(Poip)XA) = Tr{p{XA ® Ib)) (1) 

for arbitrary p Cz B and bounded operator Xa on Ha ■ If a 
state (J S H.S.{Hb) is given ouT-Ib, thenPoP = Poip)®'^ 
is a projection in B with image Bq, and is called "the 
partial trace" . Then we have a direct sum decomposition 

B = Bo®Bi, (2) 

where Bi = Pi{B) and we have put Pi = 1 — Pq. 

Now suppose we are given a hamiltonian of the form 
H ^ Ho + XHi with Hq = Ha I + 1 <E) Hb and Hj = 
X^n ^ (takc all operators hermitian). Call Z the 
generator of the one parameter group Ut of unitarities, 
on the global system B, defined by the von Ncwmann 
equation 

dtp=-i[Ho,p] = Zp. (3) 

Then we have 

Utp^ e^*p^ e-'""' pe'""*, (4) 

And one obtains the all important (superoperator) com- 
mutation property 

[Z,Po]=0. (5) 

This is telling us that the relevant degrees of freedom (the 
subsystem Bq = Po{B) do not interact with the irrele- 
vant ones (the " superpolarization space" Bi) by means 
of the unperturbed dynamics Ut- This is a fundamental 
assumption, that intrinsically links the projection Pq to 
the way the dynamics is splitted into a perturbation and 
an unperturbed one, and that should always be checked 
to avoid drastic mistakes: if this hypothesis is not veri- 
fied, the whole theory of the weak-coupling limit simply 
breaks down. 

At this point we introduce the perturbation superoper- 
ator Ap = —i[Hj,p] so that we can state with precision 
what we aim at: we simply would like to know some- 
thing about the projection at current time t of the unitary 
group describing the exact global hamiltonian dynamics. 
That is, we want to study 

pt = Po e(^+'^)'PoP = M^tV (6) 

for a given initial condition p in B. In passing, we 
have defined the subsystem evolution operator = 
Po e''^+'*"^''*Po (which of course is not a one parameter 
group of unitary transformations). 

B. Diagonal Matrices in a Closed Quantum System 

This is close to the philosophy followed by Fermi. Let 
H — C" be the n-dimensional Hilbert space over the 
complex field: then trace class operators are just n x n 
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matrices. Let these matrices form the total state space 
B, and let a hamiltonian Hq be given. Let \a) denote its 
eigenbasis, denote with Bq the space of diagonal elements 
of B in such basis, and define the projection Pq on Bq by 



(7) 



Then as before one defines Z through Hq and gives a 
perturbation A through some hermitian Hj. Obviously 
[Z, Pq] — 0, and the diagonal part of the evolution at time 
t of the initial state pinB is pt = Pq e(^+^^)*PoP = Wt^P- 



equivalently, when no there is no room for misinterpre- 
tations. 

II. TIME-SYMMETRIC MARKOVIAN 
APPROXIMATION 

Changing variables to x = s — u, a = X^u and in- 
troducing the time rescaled interaction picture dynamics 
W^-' = U_x-2rW^^2^, one is led to 



where 



1+ / U_^-2,KiX,T-a)U^-2,W^'^da, (13) 

JcT=0 



C. General Theory 

We are ready to state the general theory: we suppose 
that Pq is a (not necessarily orthogonal) projection on a 
Hilbcrt space B, put Pi ~ 1 — Pq and Bi = PiB, so that 



K{X,t)= / U-^AoiU^Aiodx. (14) 
Jo 

This form separates an "interacting" and "slowly vary- 
ing" part K{X,t) from the "rapidly oscillating" free- 
evolution to second order in the coupling constant 
A. Now in the weak-coupling limit A ^ 0, the slowly 
varying integral kernel K{X,t) converges to 



B = Bo®Bi. 



(8) 



We suppose that Z is the (skew-adjoint) generator of a 
strongly continuous one parameter group of unitarities 
Ut on B with 



UfPo = PoUt 



for all t e M, or equivalently 



[Z,Po] = 



(9) 



(10) 



and put Zi = PiZ. We suppose that ^ is a bounded 
perturbation of Z and put Aij — PiAPj. We shall al- 
ways suppose, for simplicity, that no first order term are 
present for the subsystem. That is, we shall take Aqq = 
all throughout. We let be the one parameter group 
generated by Z -I- XAn and let be the one parameter 
group generated by {Z + XA). Then putting 



(11) 



one obtains an exact and closed equation for the pro- 
jected dynamics: 



= Ut + X^ 



Ut-sAoiU^_^AiQWC dsdu. 



s=0 Ju=Q 



(12) 

This is nothing but the integrated form of the well known 
master equation constructed by Nakajima, Prigogine, Re- 
sibois, and Zwanzig 

It is a good moment here to note that, strictly speak- 
ing, Pq does not coincide with the physical subsystem: 
rather, the subsystem should be identified with the im- 
age of Pq, that is with Bq. We'll see that dissipative 
effects will be depending also upon the complete features 
of Pq itself, so not only the relevant degrees of freedom 
Bq determine the semigroup dynamics: also the way we 
look at them, with Pq, does influence the evolution. With 
this in mind, we shall use the terms "subsystem" or "Po" 



Kd= U-xAqiUxAiq dx , 



(15) 



where Kd is the celebrated Davies' superoperator. Sub- 
stituting Kd ^ K{X,t) in (|13p and moving back to the 
" Schrodinger picture" we obtain the markovian approxi- 
mation for our subsystem dynamics 



(16) 



Indeed in Q an important theorem shows that under rea- 
sonable and general conditions the approximation holds 
in the weak-coupling limit, up to A~^-rescaled times, 
independently of the subsystem dimensions or spectral 
properties. Unfortunately, although Kd is defined irre- 
spective of the subsystem spectral properties, it does not 
guarantee positivity of the generated dynamics. 

However, uniqueness of Kd, that is, of the semigroup 
approximation, is by no means to be expected: already 
in case Bq is finite dimensional (or more generally when 
Zq has discrete spectrum) one can define a time aver- 
age of Kd (which takes the diagonal part of Kd in the Z 
"super" -eigenbasis), and show that the generated dynam- 
ics is still asymptotic to the exact one in the weak-limit. 
In the infinite dimensional case this averaging map is not 
generally defined, but still no "physical" argument would 
imply uniqueness. 

With respect to this, we shall now propose a new time- 
symmetrized markov approximation: changing variable 
in eq. ([T2|) to 



(T = ^{S + U) 

r = s — u. 



(17) 



and working again in the time-rescaled interacion picture, 
we obtain 
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wy = 1 



cr=0 



r=0 



dadr. 



(18) 



r 



with g{T,(j) — |r/2— jcr — r/2||. Nowwe argue as follows: 
by direct inspection we see that in the weak-coupling 
limit A — > 0, the function X~^g(T, a), defining the inte- 
gration domain for the r-variable, converges to 



X-^9{t, a) ^ +c 



(19) 



for each tr G (0,t). But at the same time, the time- 
rescaled interaction picture dynamics moves slowly with 
respect to the r-variable (again when A '--^ 0), and this 
fact makes us guess that we can further approximate 



(20) 



in the dynamics (jl8|) . Then, the "freezed" integral in the 
r-variable factors, and our integral equation can then be 
easily seen to turn into the markovian subsystem evolu- 
tion 



with 



Kr^ I U-^AoiU^AioU-^ dx. 
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(21) 



(22) 



Here Kr is the superoperatorial and projected version of 
what Rossi proposes in [16], that also includes the second 
order energy renormalization effects. 

In fact, let us consider Zp — —i[HQ,p] and Ap — 
—i[H',p] for self-adjoint Hq and (bounded) H' on a 
Hilbert space H. Denoting A{t) = U-tAUt and H'{t) ^ 
^-iHotfji^iHot^ we factor Kr = PqKrPq where, 



Krp 



dx A 



(l)^(- 



dx 



(23) 



which in fact makes contact with the explicit form given 
in pj]. 

Clearly, starting from the interactions Aij^ the super- 
operator Kr is built in a much more time-symmetric 
fashion than K^, as can be directly seen by inspection 
of (p2| . Of course, all our arguments here are but a 
physical motivation, and we must prove that our marko- 
vian approximation indeed is consistent with the exact 
dynamics, at least up to A^^-rescaled times. From what 
we have stated here, we see that the main difficulty is to 
justify tha fact that the convergence W^fj^y^2r.^ W^'^ 

is faster than the convergence g(T, a) — > oo, when A ^ 0, 
or alternatively, that the subsystem evolves slowly com- 
pared to a A~^-time scale, in the weak-coupling limit. 



This problem presents surprising difhculties when B is 
a general Banach space, but as we have said, it can be 
attacked with a spectral analysis when we suppose B to 
be a Hilbert space, and the final answer is positive. 

The idea is precisely to compare our markov semigroup 
approximation with that of Davies (outlined before) , and 
show compatibility when A ^ 0. Compatibility of our 
semigroup with the exact projected dynamics would then 
follow automatically, thanks to the work by Davies in 
0. What we do is to estimate the greatest difference 
between the two semigroups at times from zero up to 
time A^^T, where r is some initially, arbitrarily chosen, 
positive reference time. Then we show that this difference 
goes to zero in the weak-coupling limit. So we state our 



Theorem II. 1. Suppose that 



\AoiUtAio\\ dt <oo 



Then for every r > 



lim 



sup 



o<t<x-^i 




0. 



(24) 



(25) 



We defer the proof to Appendix A, for the mathemat- 
ically inclined readers, while we now think we have com- 
pletely justified and commented the theorem from a phys- 
ical point of view. But we still have to understand the 
physical meaning of the all- important hypothesis ([24]): 
the integrand || AoiJ7t^io|| has a very simple interpreta- 
tion as the " subsystem-to-superpolarization space super- 
correlation function". In fact, reading it from right to 
left, we start in our subsystem Bo, interact and go into 
the " super-polarization space" Bi , evolve with the unper- 
turbed Ut, and finally interact a second time to land back 
into our subsystem Bq. The norm is just the mathemati- 
cal tool to evaluate how big such a "super" -amplitude is. 
Our hypothesis is thus telling us that everything works 
fine if these correlations fall faster than 1/t as the time t 
increases. That is, the super-polarization space Bi is so 
huge, that even though the global dynamics is fully co- 
herent and hamiltonian, information leaves the subsystem 
very fast, at least from a given time on. This hypothesis 
obviously excludes from a markovian approximation any 
finite (and even countably infinite) dimensional global 
system: in this case, the only things that one is going 
to see in its subsystems are nothing but some projected 
forms of never-damping Block oscillations. This explains 
the need to use an infinite-dimensional bath with contin- 
uous spectrum in Subsection II Al and the impossibility 
to observe dissipation in Subsection II BI 
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Unfortunately we are not finished yet, as the superop- 
erator Kj^ does not generate a positive semigroup. But 
the effort we have done so far in symmetrizing the well 
known usual markovian approximation Kd will be com- 
pletely paid in a moment: the idea is that Kji is, as Kd, 
highly non-diagonal in the Zq eigenbasis. But contrary 
to Kd, the upper triangular part of Kji in such a basis 
perfectly mirrors its lower triangular part. So if now we 
could somehow remove the off diagonal terms in a sim- 
ple symmetric fashion, we should be done in obtaining a 
fully diagonal superoperator. Needless to say in fact, a 
diagonal superoperator must be intrinsically linked with 
positivity features of its generated dynamics. But the 
term "diagonal part" must be handled with caution, as 
strictly speaking, it is a measure zero sector when the 
subsystem hamiltonian spectrum is continuous. We shall 
thus herepropose a "dynamic" version of the averaging 
map in that extracts a diagonal "cloud" from Kj^ 



whose width goes to zero only in the limit A 
see how. 



0. Let us 



III. SPECTRAL DIAGONALIZATION AND 
COMPLETED COLLISION TIME 

In [l| an averaging map is introduced: for an operator 
ii' : So — ^ So, we put 



= hm — 

T-.00 2T 



dq U^gKUg 



(26) 



whenever the right hand side is defined. Then Davies 
shows in 3 that if Bo is finite dimensional, then the 
operation t] is well defined and completely recovers the 
exact subsystem dynamics in the weak limit. 

Unfortunately, as said before, the spectral average t] is 
not generally defined when the system hamiltonian has 
continuous spectrum. We thus introduce a new type of 
spectral averaging: we note that the integral over time in 
eq. (p6)) can be extended to the entire real line by putting 
a T-related damping term inside the integral: 



1 



(27) 



The factor ^/^T is needed for the two definitions of [\ to 
coincide when Bo is finite dimensional: then the spectral 
decomposition 



Zo = ^ iujo 



(28) 



gives 



= lim 



dq e 



a/3 



af} 



(29) 



for both definitions ([261) and p7)). 

We may ask ourselves about the physical origin of this 
new object T: it appears that the higher its value, the 
finer is our ability to separate the system characteristic 
frequencies. This is exactly what happens when we tune 
our radio on a given frequency: the more we wait, the 
more precise we are in counting the transmitted oscilla- 
tions, and thus the finer we are in resolving the wanted 
frequency. So this T looks like a subsystem observation 
time: we have to wait an infinite amount of time to re- 
solve the system frequencies exactly, but it is everyday 
experience that we only have a finite (but long enough) 
observation time at our disposal. The idea is then to 
make use of such observation time T to compare with 
the subsystem relaxation times, and the global system 
recurrence times. It could be referred to as an "obser- 
vation time" (but also, as we shall see later on, "com- 
pleted collision time"): it is somewhat proportional to 
our instrumental uncertainty, or lack of precision, on the 
energy measurements we perform, since we study the sys- 
tem evolution only up to finite times. This way of un- 
derstanding the origin of T is consistent with the fact 
that when Zo has discrete spectrum, then the averag- 
ing map ((26)) is well defined. In fact, while we raise T, 
thus raising our energy measurements precision, we meet 
a finite value for T over which we are precise enough to 
unambiguously distinguish among the different system 
frequencies. This is not the case when the system fre- 
quencies (those of Zo) form a continuous set, as there's 
no finite value for T beyond which the system energy 
levels are distinguishable. In any case, our precision in 
measuring them gets higher as we're able to increase T . 

After having so understood the physical origin of T we 
could as well introduce a T-smoothed version of Kr (see 
ea. (P5)) for notation). 



dx 



A 



(30) 



which of course converges to the previous one in the limit 
T ^ oo, and put = PoK^Pq. Then we have yet 
another expression for the spectral average of Kn: 



K 



r> — lim Kt- 



(31) 



where we have defined 



1 

K,^-J=^j_jqe~^''-'U^,K-U, 



dtie ^Aoi(ti)/ dt2 e ^^10(^2) 



(32) 



and in the last line we have changed variable to ti — 
q + x/2 and t2 = q — x/2. This last expression deserves 
some comments: first, it is probably the most symmetric 
thing that one could imagine to build from basic objects 
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— such as Pq, Z and A — starting from the Nakajima- 
Zwanzig memory term in (I12p . Second, unUke our pre- 
hminary Kji in ((22|) , Kt gets more and more diagonal in 
the Z eigenbasis, as is evident by the procedure we use 
to obtain it. Third, the gaussian smoothing — sometimes 
referred to as homogeneous broadening — seems not to be 
there by chance, or in other words, an exponential decay 
— also referred to as inhomogeneous broadening — would 
not allow one to so easily pass from the first line to the 
second in ([5^ . In fact, the gaussian is the only distribu- 
tion such that 



Mq' + ix/2r)^mi+4)/2) 



(we shall soon be more precise on this asymptotic behav- 
ior, thanks to our next theorem). This dispenses us to 
actually compute the limit lim^^oo Kt^ that would not 
be defined in general, before performing the weak-limit 
A — > 0, and simply merges the weak-coupling and the 
collision limits together. The physical origin of T again 
helps us to understand content of this rather abstract 
procedure: when we decrease A we have to increase the 
time variable with A~^i, and so also the time T needed to 
observe the system, and resolve its frequencies, increases. 
From the mathematical point of view, constructing the 
averaging map together as the weak-limit A — > is per- 
(33) formed, allows Kt = Kt(x) to be well defined, and to 



for ti — q + x/2 and t2 — q ~ x/2. 

Now the important thing to note is that the limit 
T ^ is well defined, as said above, only when Bq is 
finite dimensional, or more generally when Zq has dis- 
crete spectrum, but, for finite values of T, Kt and Kt 
are always well defined, no matter which are the spectral 
properties of Zq nor the dimension of Bq, and moreover 
we can write Kt explicitly as: 



KtP 



1 



dtJ dt^e-^ [H'{ti),[H'{t2)]] 



= -t[H^^\p]-[£T,[^T,p]], 



(34) 



where both 



and 



2V^T 



dti dt2e~^[H'{ti),H'{t2)] 



dt 



'H'{t) 



(35) 



(36) 



are self- adjoint operators on Ti. This puts Kt in explicit 
Lindblad form, and shows that Z + }?Kt is a generator 
of a (completely positive) quantum dynamical semigroup 
(QDS) for any finite value of T (see [111). This is possi- 
ble thanks to the symmetric definition of K^: a similar 
treatment of Ko wouldn't allow to reach these conclu- 
sions, and this explains the effort we've made so far to 
justify the semigroup generated by Kr in the weak-limit. 
Of course we're not yet completely done as far as positiv- 
ity is concerned: one would like to show positivity for the 
projected Kt rather than for the unphysical Kt (we'll 
comment on this later). We shall see in the sequel that 
under fairly general conditions on Pq, Kt generates a 
completely positive QDS whenever Kt does. 

Another key issue concerns the weak-limit: if T is kept 
fixed, than of course the generated semigroup fails to 
describe the exact evolution when A ^ 0, while all of our 
approximations hold only in that case. The idea is then 
to allow a A-dependence 



T(A)-|A|~«r, A^O 



(37) 



T(A) 

generate a completely positive QDS, every step of the 
way, and as we now state, to correctly describe the exact 
dynamics in the weak-limit. 

Now we recall that we defined as 



= ^{Za+>^ KR)t 



(38) 



the semigroup dynamics generated by Kr and we give 
the following 

Theorem III.l. Suppose an "observation time" T(A) is 
given, with the asymptotic behavior reported in \31^ for 
some T > and real S^, and call 



the generated dynamics. Suppose that 

WAoiUtAwW dt < oo. 



Then for every t > 0, and < £^ < 2 we have 




sup 

o<t<A- 



(39) 



(40) 



(41) 



As with our previous theorem, we defer the proof to 
Appendix B for the mathematically inclined reader. The 
physical argument that lets us understand why we must 
suppose the asymptotics < ^ < 2 refers to some dif- 
ficulty already encountered: here we are saying that 
not only the subsystem time-rescaled, interaction picture 
evolution, moves slowly with respect to the r-variable in 
(I20p . but it also moves slowly with respect to our obser- 
vation time T(X). This obviously fixes the scale ^ < 2. 
The remaining inequality < ^ has already been dis- 
cussed when we underlined the need to make T go to 
infinity in the weak-coupling. 

As promised, we can now make some physical guess 
on the observation time T: we know that it must scale 
with |A|~^ with < ^ < 2 strictly. The most natural 
case ^ = 1 seems appropriate. Now, the meaning of T 
is that of the time needed to observe the subsystem: we 
could argue that we must observe it for a long time, if the 
interaction A is small. So we may tentatively propose 



T(A) 



1 



|A|||A|| 



(42) 
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In fact, this is to us the only "physical" time that can 
be constructed to satisfy the stated properties, starting 
from known features belonging to the original global sys- 
tem dynamics: in fact, it is clear that T must be as- 
cribed to some relaxation process, which, in our theory, 
are brought by the perturbation A — the energy brought 
by Z, instead, just makes the system (and the subsys- 
tem) rotate with unitarities — . Here the norm of the 
perturbation A has a simple meaning: it is the greatest 
eigenvalue of A. This is indeed finite, as we are supposing 
A to be bounded, and represents the greatest internal en- 
ergy transition, brought by the perturbation, that leaves 
some density matrix invariant. Then, the highest is such 
a transition, the smaller the amount of time we need to 
observe a subsystem relaxation. This explains the name 
"completed collision time": according to Eq. (|^^ . the 
more energy is brought by our perturbation, the less is 
the time required to observe the subsystem relaxation. 
This is precisely what happens with our radio: the time 
we need to reach some fixed tuning precision depends on 
how high the frequency is, that we look for. 

We have shown that the exact dynamics is well de- 
scribed up to arbitrary long times A~^r > by the 
(markovian) semigroup in the weak coupling limit. 
Such dynamics has been shown to be always well defined, 
irrespective of the subsystem spectral properties. 

IV. POSITIVITY 

We have already shown that in the (physically inter- 
esting and fairly general) case Zp = — z[_ffo,p] for a self- 
adjoint iJo on a Hilbert space H and Ap — —i[H', p] (H' 
self-adjoint), then the dynamics generated by Z+X'^Kt is 
completely positive for any given T > 0, as the generator 
itself has been written in Lindblad form: 

Ktp = ^i[H^^\p] - [Ct, [Ct^p]]. (43) 

This generator is physically completely meaningless if not 
projected, as neither it describes the global, hamiltonian, 
dynamics, nor it describes the subsystem projected dy- 
namics: we shall make use of it only for convenience in 
some intermediate step, always remembering that the fi- 
nal results only come from its projected version Kt- 

We now proceed to a simple analysis to determine 
some fairly general conditions on the projection Pq un- 
der which the superoperator L^t^ = Zq + X^Kt can also 
be put in Lindblad form, and thus generate a (com- 
pletely positive) QDS. In what follows, we shall make 
use of the dual Pq of Pq, written in Heisenberg pic- 
ture: the two objects will thus be linked by the duality 
TT{Po{p)X)=TripPo{X)). 

Theorem IV. 1. Suppose Pq is a completely positive 
projection on the observable space B{Ti.), so that it has 



Krauss decomposition 

PoX = (44) 

for some operators Va on Ti., and (possibly uncountably 
infinite) indexing setX. 

Suppose that Pq maps the observable space BiT-L) onto 
the subalgebra of observables X defined as follows: X ^ X 
if and only if for every a, 

[X, V^] = [X, V^] - . (45) 

Suppose also that Po(l) = 1, and that Pq is dual to the 
( completely positive projection ) Pq on the state space of 
trace-class operators T{H). 

Then for all real X and T > 0, the operator 1,\t = 
Zq + X^Kt generates a (completely positive and trace pre- 
serving) Quantum Dynamical Semigroup. 

Proof. We consider 

hxT=PoiZ + X''KT)Po, (46) 
and note that Z + X^Kt is dual to a generator 

LxtX = i[HxT. X] - ^{CtCt, X} + X^CtXCt (47) 

of a completely positive QDS on B{Ti), for the self- 
adjoint 

HxT = Ho + X'H^^ (48) 

(see psp and ([36|) for explicit expressions). Then for 
X E X and p G we have 

Tr{hxT{p)X) = Tr(^pPo{iZ + X''KT)X}^ 

= rr(pLAT(X)) , (49) 

for the superoperator 

Lxt = Po{Z + X^Kt) (50) 

on X. Using the fact that every X E X commutes with 
V^^^-* we compute 

Lat^ = i[J2V0HxtV0,X] (51) 

- yT.^^P^^^tVp,X} + X'J2^^^tX£tVp 

Inserting the completeness relation VdVa = 1, which 
follows from Po(l) = 1, we obtain the Lindblad form 
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a, 13 



a, 13 



.T t- 



(52) 



Q/3 



a/3 



r 



for scattering operators 



(53) 



defined on X. This shows that Lat = -^o + A^Xt gen- 
erates a completely positive QDS on So through its dual 
Lat on X: this is so indeed, precisely because 



Tr 



(e^^^* {p)X)= Tr [p e^^^* (X)) . (54) 



So this completes the proof. 



□ 



As we see, the proof is fully constructive and gives 
an explicit form for our Lindblad generator, starting di- 
rectly from fundamental objects that define the model, 
either algebraically through the projection, or dynami- 
cally through the global system hamiltonian. 

Before making some examples to show the flexibility 
of this whole theory, let us comment here on what we 
believe to be a widely spread misunderstanding: the pro- 
jected superoperator Z + ^^Kt, for example in (jiGj) . can- 
not be chosen at will. In other words, to take a general 
Lindblad operator L and project it with some given Pq, 
for studying the positivity of the dynamics e^"""*, is not 
meaningful. We have shown that it is meaningful when 
L comes from some markovian approximation that al- 
ready involves the projection Pq, through the condition 
that no subsystems first order energy renormalisation is 
present {Aqq = PqAPq = 0). It could be argued that 
this last condition is not at all necessary, and only comes 
from simplicity in the proofs involved. This is true in- 
deed, but it does not gauge the problem away: in the 
general case Ago ^ 0, the Nakajima, Prigogine, Resibois, 
and Zwanzig master equation (|12p tells us that, whatever 
the markovian approximation may be, it should have a 
projection Pq to the right as well as to the left, and an 
unavoidable Pi — I — Pq in the middle. So once again, 
this rich and intrinsic dependence between the dynamics 
and the projection cannot be causally factored: no one 
comes first. Said in other words, the (highly nontrivial) 
dynamics of interest Pqc^'^Po, is simply different from 
gPoLPo least in general. We restate that, if A^q = 
and [2', Pq] = 0, we have proven that 



is a positive map at any time if Pq fulfills the hypotheses 
of our Theorem lIV.il 

We can make a trivial counterexample along these 
lines: if one is given a global hamiltonian dynamics 
Lp — —i[Ho + XH\ p], and p is in the image of Pq, 
then PoLp = —i[Po{H), p] not only doesn't give any dis- 
sipative effect, but also (we have now abundantly shown 
that) it doesn't catch up with the exact projected dynam- 
ics Poe^^Po (that we have been studying all through- 
out), not even in the case of weak coupling limit when 
H ~ Ho -\- X^H' and A is small, and not even in case 
[Z, Po] = 0. 

This shows that the different projection techniques, 
that have been studied so far by various groups, all suf- 
fered from the dynamical restrictions present in Davies' 
theory That is, one was forced to require (a global 
hamiltonian dynamcs and) that the subsystem hamilto- 
nian has discrete spectrum. Our work here has thus been 
devoted to free oneself from this restriction. This allows 
one to study projections (i.e. subsystems) in a much more 
general contest, and we shall give an important example 
of this, straight away. We have to say that our analy- 
sis still suffers from the restriction Aqq — 0. But so did 
Davies's analysis in [1], and moreover we are currently 
obtaining already very encouraging results in the more 
general case Aqq 0. 



A. Quantum Fermi's Golden Rule 

As a first important example, suppose Va = are 
mutually orthogonal projections on the Hilbert space 7i, 
and call the "quantum populations" 



Pa = VapVa 



(56) 



for a given p E Bq- Suppose also that [Pq, Z] — 0, with 
Zp = — i[iJo, p], so that Hq — Hq, with 



(57) 



(55) 



up to A ^-rescaled times (we refer to our Theorem lIII.il 
for a more precise formulation), and the right hand side 



and put Ap = —i[H\ p\. Then the operators in ((53|) 
can be interpreted as " quantum transition rates" among 
the "quantum populations" {pa}, and one obtains an 
evolution 
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T t 



(58) 



(with i?^ ^ = V'ai/j. Vq) that in fact couples the different 
" quantum populations" /density matrices and guarantees 
positivity of each pa, as one can easily see. Since the 
projections are mutually orthogonal (VqV^ = SapVp), the 
sum over the index /3 in last equation can be restricted 
on all values such that f3 ^ a, as the condition Aqo = 
amounts to D]^^ = 0, as can easily be seen. 

The generator in ea. (|55|) for Pop, or alternatively this 
coupled linear system for density matrices pa, is of spe- 
cial interest, as it is a clear generalization of the classical 
Fokker-Planck equation, and on the other hand it is pe- 
culiar to the case Bq is not finite dimensional and Zq has 
continuous spectrum. In fact, only in these conditions 
one can hope to show that 

/•OO 

/ \\AoiUtAio\\dt<^, (59) 
Jo 

which is a necessary requirement for all the theorems on 
the weak limit that have been proven by us here and by 
Davies in P, 0]. Again, as said before, one could say 
that only when the free hamiltonian spectrum is contin- 
uous, the "polarizations" VapV/3 (a ^ P) contain enough 
degrees of freedom to allow an exponentially decaying 
solution instead of Bloch oscillations. Current work on 
the possibility of applying these ideas to the theory of 
quantum transport are being addressed by us. As it is, 
equation ([55]) could thus be addressed as a "Quantum 
Fermi's Golden Rule", as the "quantum populations" pa 
are in fact (positive definite) density matrices rather then 
(positive) real numbers. 

Note also that Eq. is a linear homogeneous system 
of (dissipative) master equations that guarantees con- 
servation of a "global trace", i.e. Tr(^^pQ,) = const, 
and positivity at all times for each quantum popula- 
tion Pa- This presents striking formal similarities with 
a very recent work Q, where the author discusses a 
most interesting "entangling" projection on bipartite sys- 
tems that could be thought of as a generalization of the 
usual partial trace projection. There, as here, the au- 
thor arrives to a linear homogeneous system of (dissipa- 
tive) master equations describing (in his case) the dis- 
sipative/decoherent interplay among the different entan- 
gled sectors that make up the subsystem density matrix. 
Now, from one side, we will show in a moment that this 
"entangling" projection can be easily incorporated in our 
theory, as it is a particular (but maybe not so much) case 
of the projections we consider in Theorem lIV.il From an- 
other side, our example here is different in that it doesn't 
refer to a bipartite system: our subsystem is, so to say, 
closed, with no environment whatsoever, and with none 



of its subsectors pa decohering to any given thermal state. 
This could implement the idea of many systems interact- 
ing together, and finally decohering to a steady state just 
because information continuously (and irreversibly) flows 
from them to their quantum polarization space. This we 
feel is much closer to the original idea that led Fermi to 
his celebrated Fermi's Golden Rule [l3|, as contrar;^to 
what has been done in recent years (see for example [Ti]), 
he did not consider bipartite systems, nor the idea that 
information irreversibly flew from a system to an envi- 
ronment. Rather, information was lost in the system po- 
larizations for some unperturbed hamiltonian basis, and 
this is exactly what we do here, extending his Fokker- 
Planck dynamics [l^ among (positive) populations to a 
Quantum Fokker-Planck dynamics among (positive defi- 
nite) density matrices. 

It could be argued that our generator is unphysical, 
or unsubstantial, as it completely disappears (it goes to 
zero) in the limit T 0. We hope to have here reported 
enough physical (and mathematical) evidence that such 
a procedure should not be undertaken (this is clear for 
example by looking at Eq. ([H])). Rather, one could, for 
example, study the steady states p = p{T) as a function 
of the collision time T, and only at that point perform the 
limit /^steady ~ limr^oo p(T). This limit will of course 
give a non-trivial result. The situation here is very sim- 
ilar to that of the weak coupling limit, at variance with 
the thermodynamical limit: the two limits simply do not 
commute, and the order to be taken depends upon one's 
interest, as justified on physical grounds. Our argument 
for the completed collision limit T — > oo is, again, that 
the semigroup approximation holds only for intermediate 
times T, between the subsystem relaxation times and the 
global Poincare recurrences (see for example or the 
brilliant exposition in (isjV Thus, it must be finite ev- 
ery step of the way, and be brought to infinity only in 
a second step. There is one notable exception to this: 
precisely when the free hamiltonian spectrum is discrete, 
there will be a finite T (as said before) over which we'll 
clearly distinguish among all the different frequencies. If 
we used a higher T' > T, the physics will then be left 
unchanged, exactly because we already reached infinite 
precision with T. This explains, physically and mathe- 
matically, the reason why one could free himself of T in 
the discrete case, by letting it go to infinity. 

So we understand why we all had to study the FGR 
in the (rather cumbersome) case of a free hamiltonian 
with mixed spectrum p^ : on one side, one needs a con- 
tinuous spectrum to hope to see dissipation, but on the 
other side, the diagonal part of the scattering operator is 
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well defined only in the discrete case. The only way out 
seems to force the interacting hamiltonian communicate 
between the discrete and the continuous part of the free 
one! Here what we do is essentially to relax the condition 
that the relevant degrees of freedom are the free hamilto- 
nian eigenvalues, and rather take them as diagonal boxes. 
Our model is thus fully preserving the quantum features 
within each sector. 



B. Entangling Projection 

As promised, we shall now make contact with the work 
in and show how the projections there considered can 
be implemented in our formalism. We just here antici- 
pate that this is going to be a generalization of the stan- 
dard trace-projection over the bath. Let H. = Ha 
be the tensor product of two Hilbert spaces, representing 
a coupled system-environment pair. Let B = T(Ti.) be 
the space of trace-class operators on H, and B{H) the 
observable algebra of bounded operators on H (the same 
goes for Ha and Hb)- 

Take a finite number C„ and Z3„ (n = 1 . . . N) , of 
operators in B{Hb), such that 

DiD„, = Snn'Bn (60) 

for (positive definite) trace-class B„ = D^^Dn in T{Hb)- 
Name the (positive definite) operators 

An = ClCn (61) 

in B{Hb) and suppose 

N 

Y^An^l, (62) 

n=l 

An An' = Snn'An , (63) 

and 

TliAnBn') ^ 5nn' ■ (64) 

Take a basis {|a)} oIHb, and define the operators 

N 

^aa' ^Ia®Y. Di\a){a'\Cn (65) 

n=l 

on H. With these operators define the completely posi- 
tive map Pq on B{7{) by 

Po^ = ^l-i,XF„„, . (66) 

OlOl' 

This map can be checked to be a projection thanks to 
the hypotheses (|60p and (|64p . The completeness hypoth- 
esis (|62p guarantees that Po(l) = I7 and hypothesis ([55]) 
guarantees that X — Pq{B{H)) is a subalgebra of the 



global observable algebra B(H). So Pq fulfills the hy- 
potheses of our Theorem lIV.il In particular, one can see 
that the projected subalgebra is 

JV 

x = {x e B{H) \x = Y,Xn®An, Xne B{Ha)} ■ 

n=l 

(67) 

Using our hypotheses we then compute 

Pop = E^-'-^^i' 

= Y,T^^B{p{l^An))<»Bn , (68) 

n 

which is precisely Equation (12) in 9]. Then it is possible 
to check that Pq is also a (completely positive) projection, 
precisely because of hypothesis . Our procedure used 
to obtain it is however different from [9j , and it would be 
surely worth the effort to compare the two approaches 
further. For example, the only common hypothesis to 
obtain such an entangling projection is our equation (j64p . 
which in [9| is Equation (6), while it seems to us that all 
other hypotheses differ in the two approaches. 

But the most fundamental difference is that here we 
have to satisfy, as amply discussed before, a dynami- 
cal compatibility hypothesis, that intrinsically links Pq 
with the form of the Lindblad operator that one wants 
to project. We have shown that if one starts from a 
free hamiltonian dynamics, Zp = —i[Ho,p] and perturbs 
with XAp = —iX[H',p], then one arrives to a Lindblad 
generator for the projected subsystem if [Z, Pq] = 0, or 
alternatively if [Ut, Pq] = for all times. For example, if, 
as normally assumed, the unperturbed hamiltonian is of 
the form Hq — Ha (8" l-f 1 (X" Hb , this dynamical commu- 
tation condition for the subsystem can be checked to be 
fulfilled if and only if for every n = 1 . . . iV we have 

[HB,An]^[HB,Bn]^0. (69) 

Then we have proven that if Aqq — our dynam- 
ics correctly describes the exact projected dynamics in 
the weak-limit, and is a (completely positive marko- 
vian) Quantum Dynamical Semigroup. If then one would 
like to take A as a Lindblad generator (instead of a 
safe self- adjoint derivation) things would get much more 
complex, as in any case one should always prove that 
— e(.Zo+>'Aoo)t jg group of isometries (see [2|), which 
gives nontrivial conditions on Aqq, when the latter is 
switched on (again a condition that links the choice of Pq 
to the dynamics). This leaves the possibility that A — Aqq 
is dissipative while Aqq is a self-adjoint derivation (i.e. it 
comes from a hamiltonian) : again, we are currently work- 
ing on this interesting and important path. But we have 
to remember that the perturbation must have a "A" at- 
tached to it, and the (dissipative) subsystem will describe 
the projected (phenomenological and dissipative) global 
dynamics, only up to second order in A, when A 0. 
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REMARKS AND CONCLUSION 



where the integral operator Ti.\ is defined by 



Although all our results suffer from being restricted to 
Hilbert spaces, we feel much confident in saying that they 
can be extended to more general (and natural) Banach 
spaces, and we are currently working on that promising 
direction. 

Another severe restriction is our hypothesis ^oo = of 
absence of first order subsystem energy renormalizations. 
Just as before, we currently obtaining very promising re- 
sults in effort to free the theory from this restriction. 

Apart from these key remarks, we have shown that, 
under the only consistency with the weak-coupling limit, 
the markovian approximation of the memory terms in the 
quantum mechanical master equation is far from unique. 
We found a generator for a Quantum Dynamical Semi- 
group that generalizes more standard approaches, and 
contrary to them guarantees both consistency in the weak 
limit and complete positivity, irrespective of the spectral 
properties of the system. We have thus been able to 
propose a new type of quantum generalization of the cel- 
ebrated Fermi's Golden Rule, giving an unprecedented 
linear homogeneous system of (dissipative) master equa- 
tions for a closed subsystem that guarantees positivity at 
all times for each of the subsystem sectors. It is worth 
noting that our model requires a continuous spectrum for 
the free hamiltonian, and thus is peculiar to our general 
theory. We have also shown that our theory dynamically 
incorporates some recently proposed models for entan- 
gled bipartite subsystems, and compared to them. All 
this opens up the way to an entirely new formalism for 
modelling nowadays mesoscopic-scale electronic and op- 
toelectronic devices, as well as to investigations on the 
fundamental, and everlasting, problem irreversibility in 
more abstract quantum theories. 



(Wa5)(t) = / dsU_x-2,KDUx-2Ms)- (A3) 
Jo 



Now H\ is a Volterra operator, so 



(A4) 



where c does not depend upon A, and we have conver- 
gence of the associated Newmann series expansion 



f^^b + Hxb + nlb 



(A5) 



So far everything can be restated for fx ^ in exactly 
the same way, leading to the corresponding definition of 
Hx. 

Subtracting the similar expansions for f and fx we 
obtain 



sup |1M^,6-W/,^6|| = ||/,-/a|| 

0<t<A-2r 



<^iiH';6-7^: 



(A6) 



It's easy to see that this series is dominated uniformly 
with respect to A, so it only remains to show that 



lim \\nx~Hx\\ =0. 



(A7) 



This means that for every e > there's a A > such that 
|A| < A imphes that for every g £V 



I dsU_x-2s^KUx-2,9{s) 
Jo 



< 



e 



(A8) 
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APPENDIX A: PROOF OF THEOREM (flTTI) 

Proof. Let V be the Banach space of norm continuous 
So-valued functions on [0,r], and let b ^ Bo- Define the 
"interaction picture" time rescaled solution 

7a (r) = U_x-^,WU,b. (Al) 
Then f^^ is a solution to the integral equation 

7x=b + nx7x, (A2) 



Before proceeding we have to keep divercences under con- 
trol: the time variable inside the integrals that define Kd 
and Kf( goes from to cxd: we show that it can be stopped 
at a finite value. To do this, we write explicitly the defi- 
nition of the operators involved in the last inequality: 



/ dsU^x-2s^KUx-2sgis) 
Jo 



dx I ds 

^0 



X-2s- 



xUx-2s- 



9is) 
(A9) 



where we have put = AqiUxAiq. Now the integral 
over the a;- variable on [0,oo] can be split into [0,x] U 
[x, cxd] , which gives 



/ ds U_x-2sAKUx^2,g{s) 
Jo 







poo 


< 




+ 2\\g\\^T \\Bx\\dx 




Jo 


J X 



13 



According to our hypothesis ([24)1 on we can choose a 
suitable x such that the last term is smaller than 
independently of A: 

2\\g\\^T j_ \\B.,\\dx<-^. (All) 

Once we fix such an our problem becomes to find a 
< A2 < Ai such that if |A| < A2, then for every < r < 
r: 



dx I ds 
Jo 



< 



We now make use of the spectral theorem: if 



iLU dE,, 



24- 
(A12) 

(A13) 



is the spectral decomposition of Z for its spectral family 
(if Zp = —i[H,p\ then the w's have meaning of en- 
ergy differences, or characteristic frequencies of the global 
system), then 



Ut 



e'"^' dE^ 



(A14) 



As for the time variable x, we split the frequency range 
into (—00, —bj) U [— tU, cJ] U (tU, 00) and define 



(TJ) 



(A15) 



Adding and subtracting U_-^_-ig_^ 
led to consider the quantities 



^'^^ and u'"^^_2g_^ we are 



dx / ds ||(C/_A-2,_, - C/i";,l.,_Ji?.C/_A-2s5(s)ll 
/O "'0 

(A16) 

and 

dx f ds\\iU_x-2,_.-U^^^_,^_^)B,U_),-2,_.g{s)\\. 

(A17) 

The first is dominated by (for compactness of the in- 
volved intervals, sup = max) 



XT max|||(C/_A-2,_^ - u'^^_2^_jB^U_x-2sg{s)\\, 
< s <T,0 < x <x} 

< XT 



+ / _ e'(^"^+^)'^ dE^ BJJ_x-2^g{s) 

-00 ■/ —UJ j 

"^(/ +/_) ii^^"^3lC^-^-'^5(s)i|. (A18) 

But since the EJs are orthogonal projections 

/oo 
\\dE^ S,C/_A-2..g(s)|| 
-00 

= \\Bjj_x-2s_g{§)\\ 

< WB.Whid (A19) 



is bounded independently of A, so we can choose ZU such 
that the quantity in the last line in (jA18|) is dominated 
by ^ independently of A (it can however depend on g, 
as a consequence of the estimation above, but this will 
not bother us): 



(^J \\dE^ B^U^x-2^g{s)\\ < ^ (A20) 



Following the same reasoning, we can choose ZU such that 
also the diflFerence in (jA17p is controlled by ^. In fact, 
again with the same reasoning, by adding and subtract- 
ing U^l^ and U^\^_j^ we easily see that we can fix an 
U such that 



dx I ds \\u[^^-,^_^B,iUx-2, ~ C^nj5(s) 
Jo 
e 



< 



27 



(A21) 



and 



dx / ds \\U[">_.^MUx-2s-^ - J.9(s 



< 



28 



(A22) 



independently of A. So we now have to prove that, when 
we fix X and tU as done before, we can find a < A2 < Ai 
such that when |A| < A2, for every < r < t: 



dx / ds 
a Jo 



< 



29' 



(A23) 

We proceed by plugging in the definition of Stiltjes inte- 
gral involved in the spectral decompositions of the uni- 
tarities Ut (see for example [l^). If S(— uJ, ZJ) is the set 
of subdivisions of the interval [— ZU, ZU], and a S E(— ZU, ZU) 
has points — ZU = cjq < t^i < • • • < ^^n„ = ZU, we put 



d{a) = max \uJk - ujk-i\ 

k—l,...,Nfj 



(A24) 



Then the limit 



N„-l 

e-' dE^f = lim ^'"""'(E^^+r - ^-J/' 

(A25) 

where ujk G [i^k,^k+i[, exists and defines the left hand 
side. Then eq. (|A23p becomes 
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lim dx ds y y 



fcl=0 fc2=0 



(A26) 



r 



with obvious notations. But clearly it suffices to consider exist a A2 > and a subdivision a such that if |A| < A2 
only one subdivision instead of two, and prove that there and d{a) < d{w) then 



/ dx ds I 



fcl.fe2=0 



29' 



(A27) 



that is, rearranging terms. 



lations by 



E 

fcl,fe2=0 



< Iwfci - Wfc2 I 



(A30) 



ds e"*^ '^('^'=l-"fc2)g(s) 



< 



29' 



(A28) 



But the left hand side is dominated by 



dx 



E 

fel,fc2=0 

- E^^^)Bx{E^^,^^, - E^^J 



(this explains the splitting in the time variable x). In 
passing, we note that we can discard the contribution 
with equal frequencies fci = ^2 in the sum for the subdi- 
vision a. 

As the rapidly oscillating integral is concerned, we take 
a A-smoothed version of the Dirac delta functional: for a 
suitable real number C > we define 



1 



(A29) 



TT X 



= ^ I dn e'"^, (A31) 



and we can safely control the slowly varying phase oscil- with fix = \M ^'^d consider the approximation 
I 



dse-'^ '"('^'=1 "'^^2)3(5) ^ / dse-'^ "Kc^fci-^fc2) / dr g{r)Sxis - r) 



2tt 



(A32) 



valid for A ^ 0. Then the term in eq. (jA29[) is bounded by 



y^ dx pfci - uJk2 1 

fci,fe2=0"^0 



2 TT 



dn 



1 



|A '^{ojkt - Wfe2) + ^\ Jo 
I 



dr WiEoj^^^, - E^^^)Bx{E^^^^^ ~ E^^Jg{r)\ 



(A33) 



Now the integral in dn can be computed easily if, for sufficiently small values for A, the integrand is never sin- 
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gular. This is the case provided that ( < 2, as we shall 
assume. Then the integrand is either always positive or 
always negative, so 



1 



Putting things together, we estimate eq. (|A29p by 

kiM=0 



-■ (A34) 



dx 



T 1 A^-C 

dx dr \\B,g{r)\\ 
h JQ 



2 TT r 



(A35) 



(the second line follows because the E^^s are orthogo- 
nal projections). Now this goes to zero as A — > 0, thus 
completing the proof. □ 



APPENDIX B: PROOF OF THEOREM (flTLT|) 

Proof. (Throughout when there's no possible misunder- 
standing we should use the notation Kt instead of the 
complete Kx(X))- With a completely analogous notation, 
we follows the very same steps in the proof of theorem 
III. II we then have to show that 



limllTiA-^All =0, 

A — *0 



(Bl) 



where both H\ and T-Cx are defined analogously to 
ea. (|A3p . respectively for and Kt{\)- Now defining 

nx by 



(Ha5)(t)= / dsU_^-2, 
Jo 

X 1^ J dq e-'i'/^'u^gKnUq^ Ux-2M, 

(B2) 



where g G V, leads to the estimation 



\nx-Hx\\ < WHx-Hx 



Hx -Hx 



(B3) 



The easiest part to estimate is the last term: one can see 
that 

\\(Mx9){T)-{'Hx9){r)\\ 



<r\\g\\ 



<r\\g\\oo dx(l-e-(t)'/^'(^))||Aoil/,^io||. (B4) 
Jo 



Now for any given x € [0, oo) one has 

lim(l-e-(t)'/^'(^)) = 



(B5) 



because of the hypothesis r(A) ~ |A|"^T with ^ > 0. 
This together with the boundedness hypothesis 



dx WAoiU^AioW < oo 



(B6) 



implies that the last term in (|B4p goes to zero as A ^ 0. 
It remains to show that 



I^A-TiAll^O, A^O, 



(B7) 



or, more explicitly, we have to show convergence to zero, 
for any g € V and uniformly on r, of 



ds U^x- 



X-^s 



X U 



A-2 



,9is)\\ 



1 ") 

dqe'^"^^"u.gKRUg-KA 

(B8) 



Now the very same arguments we used for theorem (jll.ip 
apply here, so U±q can be substituted with U±q, for a 
suitable cut-off frequency ui: the difference with eq.jBSj 
is shown to be arbitrarily small in exactly the same way 
as before. Again as before, we can evaluate the limit in 
the Stiltjes integral for the spectral decomposition over 
a single common spectral subdivision a = {—uj = loq < 
. . . < ujpf^ = oj}. So proceeding as in theorem (jll.ip we 
are now left with the problem of finding, for any given 
e > 0, a suitable A > and spectral subdivision a such 
that for any |A| < A and d{a-) < d{a) we have 



E 

kiM=0 



1 



X (i^c^fcj+i - E^^^)Kr{E^^^^^ - E^i^J 

iA"^s(i2fcj-2fc2), 



X / ds e 
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E 



1 



fel,fe2=0 

X WiE^k^+i - E^^^)Kr{E^^^^^ 

X [ ds e-'^'^'^^>'^-^>'^'^g{s) 
Jo 

< e. 



(B9) 



We multiply and divide by r|tt>fc^ — Wfcjj inside the sum 
and note that 



1 



T\uJki - ^k2 I 



1 



<C <oo 



(BIO) 



is bounded uniformly on Uk and T by a constant C. This 
means that the divergences brought by the rapid A^^ os- 
cillations cancel with the slowly varying "spectral aver- 
age" or "damping" oscillations, which we treat exactly as 



16 



in the previous theorem (see ea. (|A29p and the following, 
also for the definition of the subsequent ^) and the norm 
of the previous sum is controlled by 



T{\)nT 



N„-l 



Jo u n 



r)| 



fel,fc2=0 



A2-«C , 

< ^-l|5||_oo||i^i?,| 



(Bll) 



T(A)7r 

This shows the uniform boundedness of the last term 



with respect to \\g\\ = 1 and < t < r, so convergence 
to zero comes from the fact that 



T(A) 



A ^ 0, 



(B12) 



as we have supposed T(A) ~ |A|-«T for A ^ and ^ < 2, 
and we can take < C < 2 such that 2 + ^ — ^ > 0. This 
finishes the proof. □ 
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